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We consider theories in which there exists a nontrivial coupling between the dark matter sector 
and the sector responsible for the acceleration of the universe. Such theories can possess an adiabatic 
regime in which the quintessence field always sits at the minimum of its effective potential, which 
is set by the local dark matter density. We show that if the coupling strength is much larger than 
gravitational, then the adiabatic regime is always subject to an instability. The instability, which 
can also be thought of as a type of Jeans instability, is characterized by a negative sound speed 
squared of an effective coupled dark matter/dark energy fluid, and results in the exponential growth 
of small scale modes. We discuss the role of the instability in specific coupled CDM and Mass 
Varying Neutrino (MaVaN) models of dark energy, and clarify for these theories the regimes in 
which the instability can be evaded due to non-adiabaticity or weak coupling. 

PACS numbers: 



I. INTRODUCTION 

In order for our cosmological models to provide an ac- 
curate fit to all current observational data, it is neces- 
sary to postulate two dramatic augmentations beyond 
the minimalist assumption of baryonic matter interact- 
ing gravitationally through Einstein's equations. The 
first assumption is that there must exist cither new grav- 
itational dynamics or a new component of the cosmic 
energy budget - dark matter - that allows structure to 
form and accounts for weak Icnsing and galactic rotation 
curves. The second assumption is that a further dynam- 
ical modification or energy component - dark energy - 
exists, driving late-time cosmic acceleration. 

In the first case, recent results [l[ have added signif- 
icant support to an explanation in terms of particulate 
cold dark matter (CDM), rather than a modification of 
gravity. In the case of cosmic acceleration, however, the 
data remains consistent with a simple cosmological con- 
stant, with modifications to gravity, or with dark energy 
as the correct explanation. 

A logical possibility is that the two dark sectors - dark 
matter and dark energy - interact with each other or 
with the visible sector of the theory 0, i, i, i, i, 0, H, i, 
[13, [HI- In fact, a number of models have been proposed 
that exploit this idea to address, among other things, the 
coincidence problem [1, Q. Further, there exist classes of 
modified gravity models which may either be mapped to 
interacting dark energy models, or closely approximated 
by them over a broad range of dynamical interest [ill , fl^ . 

There are concerns, however, about coupling these two 
rather differently behaving sectors. One concern is the 
presence of dynamical attractors that produce a cosmic 
expansion history significantly different from ACDM; 
see, for example, Refs. [l,[i3,[lJ]- Another specific exam- 
ple is the possibility of instabilities that are not present 
for the uncoupled system [l^, [l^, [13 . 

In this paper we perform a careful analysis of the vi- 



ability of coupled dark energy-dark matter models. We 
consider implications of such coupled theories on cosmo- 
logical scales, both in terms of homogeneous background 
expansion and the growth of linear perturbations. We 
focus, in particular, on a class of models in which there 
exists an adiabatic regime in which the dark energy field 
instantaneously tracks the minimum of its effective po- 
tential, as explored by Das, Corasaniti and Khoury[8[. 
We show that if the coupling strength is much larger 
than gravitational, then the adiabatic regime is always 
subject to an instability. The instability is characterized 
by a negative sound speed squared of an effective cou- 
pled dark matter/dark energy fluid, and results in the 
exponential growth of small scale modes. We analyze 
a number of different models, and show that for these 
models the instability strongly constrains the region in 
parameter space that is compatible with observations. A 
short version of our results was given in the recent paper 

We imp rove on previous investigations of this insta- 
bility [l5l . [TgI . [r^ in a number of ways. First, we show 
that the instability occurs only for coupling strengths 
that are strong compared to gravitational coupling, and 
can be evaded at weaker couplings; this point was missed 
all in previous work. Second, we give a simple intuitive 
explanation of the instability as a type of Jeans insta- 
bility. Normally, for cosmological perturbations, Hubble 
damping converts the exponential growth of gravitation- 
ally unstable modes into power law growth. However, 
here the Hubble damping is ineffective, due to the fact 
that the effective Newton's constant for the interaction of 
dark matter with itself is much larger than the Newton's 
constant governing the background cosmology. The re- 
sult is the exponential growth of perturbations. Finally, 
we generalize previous treatments to allow an arbitrary 
coupling between the dark energy sector and the visible 
sector, in addition to the coupling between dark energy 
and dark matter. 
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In more detail, we consider models which are charac- 
terized by a function ac{4>) governing the interaction be- 
tween CDM and a quintessence type scalar field cj), and a 
function ab{4>) governing the interaction between visible 
matter (baryons) and cf). The effective Newton's con- 
stants for the interaction of CDM with itself (cc), the 
interaction of CDM with baryons (cb), and the interac- 
tion of baryons with themselves (bb) are (see Appendix 



Gcc - G[l + 2mla'M^] 

Gcb - G[l + 2mla'Mab 
Gbb - G[l + 2mla',{cbf] 



(1.1a) 
(1.1b) 
(1.1c) 



Here rrip is the 



These Newton constants are those of the Einstein frame 
in the short wavelength limit, and include the effect of 
the scalar interaction mediated by 
Planck mass. 

The adiabatic instability occurs when Gcc S> G, or 
equivalently rupa'^ 3> 1, when the scalar coupling be- 
tween dark matter particles is strong compared to the 
tensor coupling^. In the present day Universe, this 
regime is excluded by observations if we assume that per- 
turbations to the cosmological background value of (j) are 
in the linear regime, so that the parameters (jl.ip are 
constants (see below). First, tests of general relativity in 
the Solar System constrain mpaj, to be small compared to 
unity. Second, observations of tidal disruption of satellite 
galaxies of the Milky Way provide the constraint^ jTo| 



Gbc 



\J GccGbb 



- 1 



< 0.02. 



(1.2) 



Combining these constraints excludes the regime Gcc 3> 
G. 

Nevertheless, it is still of interest to explore the adia- 
batic instability, for a number of reasons. First, the ar- 
gument given above assumes a single type of CDM parti- 
cle. In more complicated models with two different CDM 
species, one coupled and one uncoupled, the constraint 
can be evaded (see Sees. IV B] and IVDl below). Second, 
the Newton's constants can evolve as a function of red- 
shift, via the dependence of </> on redshift, so present day 
observations (z w 0) do not exclude the occurrence of 
the instability at high redshifts. Third, it is useful to 
have multiple independent constraints on coupled mod- 
els. Finally, it is not necessarily true that perturbations 
to the cosmological background value of (f> are in the lin- 
ear regime. In particular, it has been claimed that obser- 
vational constraints on mpa^ can be evaded in chameleon 



models^, due the dependence of the local effective New- 
ton's constants (|l.ip on cj) and hence on the local matter 
density 

The structure of the paper is as follows. In the next 
section, we introduce the general class of models and 
the framework in which we work. In section IIIII we dis- 
cuss the adiabatic regime, giving both local and nonlocal 
conditions for its applicability. We also derive the effec- 
tive equation of state for these models in the adiabatic 
regime; as previously noted in Ref. Q, superacceleration 
is a generic feature of these models. Section HVl discusses 
the adiabatic instability that can arise in the adiabatic 
regime. We discuss two complementary ways of under- 
standing the instability, a hydrodynamic viewpoint and 
a Jeans instability viewpoint. We also derive the range 
of lengthscales over which the instability operates, which 
was incomplete in earlier analyses [l^, . In section |V] 
we apply our general analysis to some well-known coupled 
models. We perform analytic and numerical analyses of 
the evolution of the background cosmology and of per- 
turbations in both coupled CDM, chameleon [l^,[23| and 
mass varying neutrino (MaVaN) t22„^, models, and 
identify the regimes in which they are subject to the in- 
stability. We conclude in I VII bv summarizing our findings 
and discussing their implications. 

In appendix |^ we generalize the derivation of the in- 
stability to include baryonic matter, and in appendix [Bl 
we generalize the derivation of the instability to regimes 
where CDM cannot be treated as a pressureless fluid but 
instead must be treated in terms of kinetic theory. Fi- 
nally, in appendix [C] we derive and describe the Jeans 
instability viewpoint on the instability. 

A note on conventions: throughout this paper we use 
a metric signature (-,+,+,+), and we define the reduced 
Planck mass by TOp = (SttG)" 
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II. MODELS WITH A COUPLING BETWEEN 
DARK ENERGY AND DARK MATTER 

A. General class of models 



We begin from the following action 

S = 5[.9a6, 



.[p2aj(0) 



9ab,'^j], 



(2.1) 



This requirement is admittedly a fine tuning and unnatural from 
the viewpoint of effective field theory. 

Only the combination 111.21 1 of the Newton's constants can be 
constrained by observations of the gravitational interactions of 
baryons and dark matter, since once cannot separately measure 
the densities and the Newton's constants. 



where gab is the Einstein frame metric, is a scalar field 
which acts as dark energy, and ^I'j are the matter fields. 



In particular, chameleon models in the regime m.pa'^ 2> 1 have 
been extensively explored [21II . 
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The functions aj{(j)) are coupling functions that deter- 
mine the strength of the coupling of the jth matter sector 
to the scalar field. The matter sectors include cold dark 
matter with coupling function Q!c(0), and baryons with 
coupling function ab(0). 

This general action encapsulates many models stud- 
ied in the literature. The case of equal coupling to dark 
and visible matter, Q;j(0) = a{(/)) for all j^, corresponds 
to scalar-tensor theories of gravity which have been ex- 
tensively studied, recently under the name of "coupled 
quintessence" This class of theories includes the /(i?) 
modified gravity theories 25[ . Several authors have con- 
sidered the case ab{4>) = 0, in which the dark energy 
couples only to dark matter. We note that this choice 
yields a microphysical model for Modificd-Source Grav- 
ity [2^ in the adiabatic regime discussed in Scc. lIIII bclow 
(in the approximation where one considers only a dark 
matter source and neglects baryons). 

There is good theoretical motivation for considering 
nontrivial and different coupling functions Q:j(0), since 
this is a generic prediction of string theory and of higher 
dimensional models. In fact, typically the moduli and 
dilation fields of string theory must be massive today, 
because for massless fields it is difficult to satisfy the ob- 
servational constraints [Solar System tests and fifth force 
experiments that require dab/d<f> < 10~'^m~^ , equiva- 
lence principle tests that require matter coupling func- 
tions for different matter sectors aside from dark matter 
to be very nearly the same] in a natural way because of 
loop corrections. However dynamical dark energy models 
require massless or nearly massless fields so one is forced 
to confront the naturalness issue. One possible solution, 
suggested by Damour and Polyakov [23|, is that there 
is an attractor mechanism under cosmological evolution 
that drives the theory to be very close to general relativ- 
ity, aj for all j. It is possible that this mechanism 
does not work perfectly, and that there arc residual de- 
viations in the form of nontrivial matter couplings. 

In addition, while equivalence principle tests strongly 
constrain differences between the coupling functions for 
different types of visible matter, the corresponding con- 
straints on dark matter are much weaker, as first pointed 
out by Damour, Gibbons and Gundlach So there is 
good motivation for exploring models with various cou- 
plings to dark matter. 

If the jth sector consists of a fermion of mass m, the 
corresponding action in Eq. (j2.ip can be written as 



been explored in the literature. For example, Farrar and 
Peebles [l^ consider the case p = 0, and Bean consid- 
ers the case p ~ q. All of these choices are equivalent^ in 
the nonrelativistic limit, when the fermion's rest masses 
dominate their gravitational interactions. The choice we 
make here is motivated by the fact that it satisfies the 
equivalence principle within the jth sector, and also the 
fact that it arises naturally from higher dimensional mod- 
els. 

The field equations resulting from the action (|2.ip are 

mlGab = Va0V60 - iga6(V0)2 ~ y(0)5^b 

+ J2 e'"^^^^ [(« + aU^ b + mab] , (2.3) 



and 



V,V>- y'(0) = 5]a;(</>)e4"^(^)(pj - 3pj). (2.4) 



where the prime represents a derivative with respect to 
(j). Here we treat the matter as a fluid with Jordan- 
frame density pj and pressure pj and with a 4- velocity Uj a 
normalized according to g'^^u-^aU^b = —1 [see Appendix 

For much of this paper we shall work in the approxi- 
mation where we consider only the gravitational effects 
of the dark matter and neglect the gravitational effects 
of the baryons. In this approximation, the sums over j in 
Eqs. (|2.3p and (|2.4p reduce to a single term involving the 
dark matter. Dropping the index j for simplicity, and de- 
noting the CDM coupling function ac{4>) simply by a{4>), 
the field equations become 

mlGab = Va?iVfc?i - iga6(V<?!))2 - y(0)5^,, 

+e4"(^)[(p + p)ua«f,+raafc], (2.5) 



and 



VaV^^ - F'(0) = a'((/))e4"(^) {p - 3p), (2.6) 



These are are the standard equations for a scalar-tensor 
cosmology. 

For dark matter we have p = 0, and we define a 
rescaled density variable 



qa( 



-'mvf'iVE'i 



(2.2) 

with p = 2 and q = 3. Other possibilities for p and q have 



_ 3a(d>) - 

p = e p. 



With this new variable, we obtain 



1 



(2.7) 



m;Gab = Va0Vb0 - -gabi'^ckY - V{(t))gab + e°'^f>pUaUb, 

(2.8) 



A more precise characterization of equal coupling is a'j{4>) = 
a' {(f)) for all j, since any constant term in aj{<f)) can be absorbed 
into a rescaling of all dimensionful parameters in the action Sj . 



Up to redefining a by multiplying by a constant. 
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and 

V„W - V'{cj>) = a'{^)e^^^'^p . (2.9) 

The scalar field equation (|2.9p can also be written as 

V,V>-Kff('/') =0, (2.10) 

where the effective potential that includes the matter 
coupling is 



l/eff(0) = n'^)+e"('>V- 



(2.11) 



and the prime means derivative with respect to (f> at fixed 
p. The equations of motion for the fluid are [see Ap- 
pendix [K\ 



and 



Vaipu'') = 0, 



(2.12) 



(2.13) 



Note that the equations of motion (|2.8p - (|2.13p of 
the theory do not depend on the coupling function ai,{4') 
of the scalar field to visible matter, in the approxima- 
tion where we neglect the baryons. Our discussion of the 
instability in the following sections will be valid for arbi- 
trary ab{4>)- However we note that ab{4>) will enter when 
we want to compare with observations, since the metric 
that is measured is e'^°"'gab- When discussing observa- 
tions below we will focus on the case ab ~ 0, where the 
dark energy is coupled only to dark matter. 



III. THE ADIABATIC REGIME 



The effective potential (|2.1ip governing the evolution 
of the scalar field is the sum of two terms, one arising from 
the original potential T^(</>), and the other arising from 
the coupling to the energy density of the dark matter 
fiuid. It is possible, for appropriate choices of the poten- 
tial and the coupling function, that competition between 
these terms leads to a minimum of the effective poten- 
tial. Further, in some regimes, it may be possible for 
the solution of the equation of motion for </> to adiabati- 
cally^ in the track the position of this minimum. That is 
to say, the timescale or lengthscale for cj) to adjust itself 
to the changing position of the minimum of the effec- 
tive potential may be short compared to the timescale or 
lengthscale over which the background density is chang- 
ing. 

This adiabatic regime has been previously discussed for 
spatial variations of (f> in the interior of massive bodies 



Note that throughout this paper, we use "adiabatic" in the sense 
of "gradual change" , not in the thermodynamic sense of "at con- 
stant entropy". 



in the so-called Chameleon field models [1^, [20| , and for 
the time variation of (f) in a cosmological context by Refs. 
[1, [lB|. It also has been studied for the specific case 
of f{R) modified gravity models [2^, [s^, [3l|, for which 
the action in the Einstein frame is of the general form 
(12. ip . In this section we review the adiabatic regime in a 
general context and give a careful derivation of its domain 
of validity. 

The adiabatic approximation consists of (i) omitting 
the d'Alembertian term in Eq. (|2.10p . which gives an al- 
gebraic equation for </> that one can solve to obtain ^ as a 
function of the density p; (ii) omitting the terms involv- 
ing the gradient of </> from the field equation (|2.8p . The 
resulting equations are the same as those of Modificd- 
Sourcc gravity p6j , in the approximation where one con- 
siders only a dark matter source and neglects baryons. 
The equations can be written as general relativity cou- 
pled to a fluid with an effective energy density peft and 
effective pressure Pctf- 



rripGab 



where 



PcS 



[{PoS + Pes)UaUb + PeSgab] , (3.1) 



(p) = e"[*-('')V + ^['/'m(p)] , (3.2) 



pM = -V[<l>m{p)] ■ (3.3) 
Here 4>m{p) is the solution of the algebraic equation 

K'ff(0) = y'(0) + a'(<^)e"(«p = O (3.4) 



for (f). Eliminating p between Eqs. p.2p and (|3.3p gives 
the equation of state Peff = Pcff (poff). 

In the adiabatic regime, the matter and scalar field are 
tightly coupled together and evolve as one effective fluid. 
By taking the divergence of the flcld equation p.ip . we 
see that this fluid obeys the usual fluid equations of mo- 
tion with the given effective equation of state. In a cos- 
mological context, the effective fluid description (|3.ip is 
valid for the background cosmology and for linear (and, 
indeed, nonlinear) perturbations. Therefore, the equa- 
tion of state of perturbations is the same as that of the 
background cosmology, or the so-called entropy pertur- 
bation vanishes. 



A. Condition for global validity of adiabatic 
approximation 

Of course, the adiabatic approximation will not be a 
good one for all choices of V{4>) and a{(j)), or in all phys- 
ical situations. Our goal in this subsection is to establish 
criteria under which we can trust the adiabatic approx- 
imation. Roughly speaking, the mass of the scalar field 
associated with the minimum of the effective potential 
must be sufficiently large. More precisely, we define the 
effective mass as a function of density p by 



Ap) 



(3.5) 



= 0m(p) 
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where the derivatives are taken at constant p, and (t>mip) 
is the value of the scalar field which minimizes the effec- 
tive potential at a given density p. We assume that mlg 
is positive, otherwise there is no local minimum of the 
potential and an adiabatic regime does not arise. For a 
perturbation with lengthscale or timescale C and density 
p to be in the adiabatic regime, it is necessary that 



We call this condition the local adiabatic condition^ If 
the condition is satisfied everywhere in spacetime, then 
the adiabatic approximation will be good everywhere. 
However if the condition is satisfied in a local region, 
it does not necessarily follow that the adiabatic approx- 
imation is valid in that region, due to non-local effects 
ESIlll. This is discussed further in Sec. [mB] below. 

We now derive a slightly more precise version of the 
local adiabatic criterion p.6p . For a given density po we 
define 0o = (^m(/3o) and 5(f) = — ^o- We expand the 
potential as 



V{(t)) = + ViSd) + ^^2^0' + 0{5f) , 
and defining W{(j)) = e"*^"^^ we similarly expand 

W{(j)) ^Wo + WiScf) + ^W254)^ + 0{S(t)^) . (3.8) 
The effective mass ()3.5p is then given by 

i)Q,Po) ^V2+ paW2 , (3.9) 



'eff 



and the condition (|3.4p that the effective potential be 
minimized yields 



(3.10) 



= -poWi. 



Also a short computation using the definition (|3.4|) of the 
function (j)n\{p) gives 



dp 



"T-off(Po)^ 



(3.11) 



where ao = a(?l'o) and aj, = a'((/)o). We also define L to 
be the smallest lengthscale or timescale over which the 
density p changes, so that \l a^"'P ^ pj and (Vp)^ ~ 
P^IL\ 

The field equations arc the Einstein equation (|2.8p . 
the scalar field equation (|2.9p . and the fluid equations 
(|2.12p and (|2.13p . These equations are not all indepen- 
dent, since the Einstein equation enforces conservation 
of the total stress energy tensor. We will take as the in- 
dependent equations just the Einstein equation (j2.8p and 



the fluid equations (|2.12p and (|2.13p . since the scalar field 
equation can be derived from these. ^ 

Therefore, to justify the adiabatic approximation, it 
is sufficient to justify dropping the scalar field derivative 
terms in the Einstein equation (j2.8p . The ratio of the 
scalar field gradient terms to the potential term evaluated 
at the adiabatic solution is of order 



(3.6) {y^f 



1 



(Vp)^ 



1 



V{(i>) Vo\dpJ ''''' Vt:>L'^\d\up 
By combining Eqs. (|3.10p and p. lip we obtain 

dlnp m^jj' 



(3.12) 



(3.13) 



and using this to eliminate one of the factors of In p 

from Eq. ((37T2)) gives 



.(3.14) 



(V0)2 _ Vl d<p^ 




d\nV 




y(</>) yodlnp' 




dlnp 





(3.7) Thus, the adiabatic approximation will be valid whenever 



1 



d\nV 



d\np \TO^jj£2 



<C 1. 



(3.15) 



Now, since the first factor on the left hand side involves 
derivatives of logarithmic factors, we expect this prefac- 
tor to gencrically be of order unity. When this is true the 
condition (|3.15p reduces to the local adiabatic condition 
(EH). 



B. Nonlocal condition for breakdown of adiabatic 
approximation 

As mentioned above, the condition L 3> rii^^ in a lo- 
cal region is a necessary but not a sufficient condition 
for the validity of the adiabatic approximation in that 
region. This is because the corrections to the adiabatic 
approximation are determined by a wave equation ob- 
tained by perturbing Eq. (|2.9p whose solutions depend 
in a non-local way on its sources. The corrections in a 
given region can become large due to a breakdown of the 
local adiabatic condition (|3.6p that occurs elsewhere [3l| . 

For the special case of static, spherically symmetric 
systems, and for a constant density object in a constant 
density background, the chameleon field papers [TH, 
derived a precise condition for the validity of the adi- 
abatic approximation inside the system, the so-called 
"thin shell" condition, which depends on the asymptotic 
value of the potential. The thin-shell condition is both 



It has also been called the Compton condition [3lj , since the RHS 
of Eq. 113. 6| I is the effective Compton wavelength of the field. 



Similarly, the approximate form 113.41 1 of the scalar field equation 
can be obtained from the approximate form l|3.1ll of the Einstein 
equation together with the fluid equations. 
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a necessary and sufficient condition, but it restricted to 
static situations. Below we will show that for several 
specific examples the thin-shell condition and the adia- 
batic condition (|3.6p give the same predictions in order 
of magnitude for the boundary of the adiabatic regime. 

Wc now discuss an order-of-magnitudc non-local cri- 
terion for static, spherically symmetric situations that 
predicts when the adiabatic approximation breaks down 
even when the local adiabatic condition (|3.6p is satisfied. 
The condition is a generalization of a condition derived 
by Hu [3l| in the context of f{R) modified gravity mod- 
els, and is also a generalization of the thin-shell condition 

of Refs. [ii,!!^. 

Consider a spherically symmetric density profile p(r), 
which we will assume for simplicity is monotonically de- 
creasing as r increases. From this density profile wc can 
compute the corresponding adiabatic scalar field profile 



(f>i,dir) = (f>m[p{r)]. 



(3.16) 



If the local adiabatic condition (|3.6|) is satisfied for all r, 
then this adiabatic field is a good approximation to the 
actual solution (/)(r). Suppose therefore that the local 
adiabatic condition is violated^ for some set of values 
fi < r < f2. Then typically what occurs is that there is 
a larger interval ri < r < r2 with ri < fi and r2 > r2 
on which the solution 0(r) differs significantly from the 
adiabatic field profile p.l6p . For a given interval (ri, r2), 
significant deviations from the adiabatic approximation 
should occur in the vicinity of r = ri if 



> 



n 



[p{r) - P2] , 



(3.17) 



where pi = p(n), = (jJadiri), ai = a{(t)i) and a- = 
a' {(pi), i — 1,2. The criterion assumes that 7*2 3> ri and 

P2 < Pi- 

The criterion (|3.17p is a generalization of other criteria 
that have appeared in the literature. First, if one assumes 
that the density varies on a Icngthscale of order ~ r, then 
using this to approximate the integral in (j3.17p gives 



> 



e"'rl{p,~ p2). 



(3.18) 



Equation (|3.18p is the condition derived in Rcf. [31I in 
the context of f{R) models, although there the factor of 
a'^e"'^ was neglected. The condition p.lSp was found to 
reliably predict the onset of deviations from the adiabatic 
profile in numerical solutions for the Solar System and 
the Galaxy 

Second, one can derive from Eq. (|3.17p the thin-shell 
condition of Refs. [l^, [13], up to a factor of order unity. 
Suppose one has a uniform density sphere of density pc 
and radius Rc, embedded in a uniform density, infinite 



® If the density goes to a constant at large r, then the local adia- 
batic condition is always satisfied as r — > 00. 



medium of density Poo- For this case the local adiabatic 
criterion is satisfied everywhere except at the point of 
discontinuity of the density at r = Rc, so we expect 
a breakdown of the adiabatic approximation near this 
point. Following Ref. (l9| we assume that a' = /3/mp is a 
constant, and that we are in the regime where q;(0) ^ 1, 
so e" 1. If we apply the criterion (|3.17p to an interval 
(ri,r2) with ri < Rc < r2, wc obtain that the adiabatic 
approximation should fail near r = ri if 



A(t) 



> 



i?3 



riRl 



(3.19) 



where $c ^ PcRl/ f^p is the Newtonian potential at the 
center of the sphere. If the left hand side is large com- 
pared to unity, then Eq. (|3.19p will be satisfied for all 
values of ri except very close to = 0, and it follows 
that the adiabatic approximation will not apply in most 
of the interior of the sphere. On the other hand, if the 
left hand side of Eq. (|3.19p is small compared to unity, 
then it follows that one would expect the adiabatic ap- 
proximation to be valid throughout the interior of the 
sphere except very close to the surface, in a thin shell of 
thickness ~ RcA(f>/{f3mp^c)- Both of these conclusions 

agree in order of magnitude with those of Refs. [H, H^l 
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We now turn to the derivation of the non-local condi- 
tion (|3.17p . We start by noting that static solutions of 
the equation of motion (j2.10p can be obtained by extrcm- 
izing the energy functional 



E= dr 



p{r) 



(3.20) 



The basic idea [3l| is that the adiabatic field profile 
p.l6p minimizes just the potential energy, and there 
some kinetic energy cost for following this adiabatic pro- 
file. When this kinetic energy cost becomes sufficiently 
large, it becomes energetically favorable for the field to 
switch to a different, non-adiabatic profile with a smaller 
kinetic energy and with a larger potential energy, for a 
net gain in energy. 

We now compute the total energies Es.d for the adi- 
abatic profile and i^triai for an alternative trial profile 
which is qualitatively similar to numerical solutions [sij . 
When 



-Atrial — £^ad < 



(3.21) 



there is a net gain in energy for switching to the trial pro- 
file, and so the adiabatic profile is no longer a good ap- 
proximation. The trial profile is simply Atrial ('') = ^2 = 
constant for ri + Ar < r < '"2, and Atrial (?") = a — (3/r for 



This example shows that the condition 113.181 derived in Ref. [3111 
is less general than the condition Il3.17( derived here, since that 
condition cannot be used to derive the thickness of the thin shell. 
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ri < r < r+Ar, with a and f3 chosen to satisfy continuity 
at ri and ri+Ar, </)triai('^i) = (^i and (/)triai(fi + Ar) = (/)2. 
This trial profile depends on one parameter, namely the 
width Ar of the region in which the field transitions from 
01 to (j>2- We obtain 



trial 



ad 



1 



0ad(O' 



+a^e"^ [0triai(r)-0adM] [p(r) 



ri(ri + Ar)(02 -(^i)^ 



2Ar 



P2j 



(3.22) 



In deriving this formula we have used an approximate 
Taylor expansion of the effective potential about the 
point 02 together with Eq. (|3.4p : 



90 



1,' / 



-0)(P-P2).(3.23) 

Consider first the kinetic energy contributions, the first 
and third lines of Eq. p.22p . The minimum kinetic en- 
ergy for any field configuration which satisfies 0(ri) — 0i 
and 0(r2) = 02 is the energy of the Laplace equa- 
tion solution with these boundary conditions'' , namely 
£;K.mm = (A0)Vir2/(2(r2 - ri)), where A0 = 02 - 0i. 
Therefore we can write the kinetic energy of the adia- 
batic field profile [the negative of the first line of Eq. 
(|3.22[) ] as i^K.ad = £^K,min(l + £) , whcrc e is dimension- 
less and nonnegative. We will assume that e > 1, since 
e ^ 1 would require the adiabatic profile 0ad(f) to be 
very close to the 1/r profile, which would be a fine tun- 
ing. For generic density profiles we expect e > 1. If we 
now choose Ar — 2ri/e, then it follows that the net gain 
in kinetic energy is 

^E]S^>E^^^,^^r,Aq^\ (3.24) 

using r2 > ri . 

Turn now to the potential energy contributions, the 
second line of Eq. (|3.22p . If we assume that the potential 
V{(f>) is monotonically decreasing, as in the Chameleon 
field models, then it follows that 0ad(?') is an increasing 
function of r, so that 

0trial(7') - 0ad(?') < 02 - 01- (3.25) 

This implies that the potential energy term in Eq. (|3.22[) 
is bounded above by 



A0 / [p{r) - p2] 



Inserting this together with the estimate (|3.24p of the 
kinetic energy gain into Eq. (j3.2ip yields the criterion 



C. Observed equation of state parameter 

In this section we derive the observed equation of state 
parameter iWobs for these models, in the adiabatic regime, 
for the case a;, = 0. Generically we have Wohs < — 1 
corresponding to superacceleration, as previously noted 
in Ref. Q. 

The background cosmological evolution is given in the 
Einstein frame by the equation 



(3.26) 



where po is a constant, together with the evolution equa- 
tion for the scalar field. Observations of the acceleration 
of the Universe are fit to the model 



PDE(a) 



(3.27) 



where pi = e"°po is the observed matter density today'^, 
aQ is the value of a today, and pde is the inferred "dark 
energy density" . The equation of state parameter is then 
given by 



Wobs(a) 



-1 



1 dlnpDE 
3 din a 



(3.28) 



Combining Eqs. (P?^ ~ ((5:^ and using Eq. ([23) to 
eliminate the term proportional to d(j)/da gives 



■U^obs 



-V 



(3.29) 



Next we use Eq. (|3.4p again to obtain po/a? 
— e~"T^'(0)/a'(0), which gives 



Wohs 



1 



din V 



i^(l 



(3.30) 



This formula has the property that Wohs = — 1 today 
for all models. Also expanding to first order about a = 
1 we obtain l/wobs = — 1 + ln(y/Vb), where Vb is the 
value of V today, so i«obs < — 1 in the past since V — 
F[0m(/3)] > Vo in the past. 



IV. ADIABATIC INSTABILITY 

In the adiabatic regime, the models (|2.ip discussed 
here can exhibit instabilities on small scales character- 
ized by a negative sound speed squared of the effective 
coupled fluid. This instability extends down to the small- 
est scales for which the adiabatic approximation is valid. 



Note that in this model some fraction of the mass density in the V 
term in Eq. 113.261 1 will cluster, so measurements of mass density 
This can be shown using the Schwarz inequality (/ fgdr)^ < using clusters will not correspond exactly to measurements of the 

/ f^dr J g^dr, with / = r<f>'{r) and g = 1/r. second term in Eq. 113.261 1. We will not address this issue here. 
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Starting with a uniform fluid, the instabihty will give rise 
to exponential growth of small perturbations. The final 
state of the coupled fluid is beyond the scope of this pa- 
per. Theories that exhibit this instability are typically 
ruled out as models of dark energy. 



that the flrst derivative of VcS with respect to 4> vanishes 
at (j) = 4)^(p), we obtain 



da 



dcP' ' 



(4.4) 



A. Hydrodynamic viewpoint 



where the derivatives are taken at constant p. Using Eq. 
(|2.1ip and then eliminating p using Eq. (|3.4p gives 



This adiabatic instability was first discovered by Af- 
sliordi, Zaldarriaga and Kohri in a context slightly 
different to that considered here. That context was the 
mass varying neutrino model of dark energy, where a dy- 
namical dark energy model is obtained by coupling a light 
scalar field to neutrinos but not to dark matter. The in- 
stability was previously discussed in the context of the 
models considered here by Kaplinghat and Rajaraman 
[l^. In this paper we will generalize the treatment of 
the instability given in Ref . [16| . 

A fairly simple form of the stability criterion can be ob- 
tained by writing the potential V^(0) as a function V{pi) 
of the coupling function a((/)) by eliminating 0. This gives 



Pcff = + e"p = 1/ - = ^ - 3- ' 

daldq) da 



(4.1) 



where we have used Eqs. p.2[) and (|3.4p . The square 
of the adiabatic sound speed is then given by 



1 _ dpcff _ dp^s/da _ £ [V" - ^] _ ^ ^ 
dpcs dpcs/da ^ [-V] 



dV 
da 



The system will be unstable if 



cl <0. 



(4.2) 
(4.3) 



Equation (j4.2p gives a simple prescription for computing 
when a given theory will be unstable, assuming it is in 
the adiabatic regime. This equation furnishes as a 
function of a, which can be re-expressed as a function of 
(j) using a = a{(t)), and then as a function of density p 
using = (j)ra{p) from Eq. (|3.4p . 

More generally, outside of the adiabatic regime, the 
effective sound speed is the ratio of the local pressure 
and density perturbations 

'~ 5p 

and can differ from c^. In a cosmological context we 
have = P/p and Cs(fc,a) = 5P{k,a)/5p{k,a), where k 
is spatial wavenumbcr and a is scale factor. For most of 
this paper we will consider only the adiabatic regime in 
which — > c^, although the more general regime will be 
probed in our numerical integrations in Sec. |Vl 

We now argue that the adiabatic sound speed (|4.2p is 
always negative in the adiabatic regime. In the definition 
(|3.5p of the effective mass, we rewrite the (p derivative in 
terms of a derivatives using a = a{(j)). Using the fact 





-d^v 


dV 




da"^ 


da 



(4.5) 



Simplifying using the expression (|4.2[) for the sound speed 
squared gives 



da 



dV 1 

da ci 



(4.6) 



The adiabatic instability arises in the case < and 
TO^ff > 0. [When m^g < there is no adiabatic regime, in 
the sense that adiabatic solutions arc subject to a tachy- 
onic instability. This instability has been discussed in 
the context of f{R) gravity models in Refs. [s^. l33l. [sij .] 
It follows from Eqs. (|4.2p and (|4.6p that the adiabatic 
instability occurs when 



or when 



— < — 

da da^ 



— < — 

da da^ 



(4.7) 



(4.8) 



However, the region in the (T^q, Kqo) parameter space 
defined by the union of the regions (|4.7p and (|4.8p is 
just the adiabatic regime; the adiabatic approximation 
requires both that V^a < so that the effective poten- 
tial has a local cxtremum, and also that m^j^ > so 
that the local extremum is a local minimum. It follows 
that the adiabatic instability condition (|4.3p is satisfied 
throughout the adiabatic regime. This gencricity of the 
instability has also been deduced by Kaplinghat and Ra- 
jaraman [Toj using a different method. However, as we 
discuss in the following subsection, the instability actu- 
ally occurs only if the coupling oi'{(f)) is sufficiently large, 
a point missed in Refs. [la, [iF 



B. Scales over which the instability operates 

For the instability to be relevant on some spatial scale 
C, then we must have 



C > mj 



(4.9) 



in order that the adiabatic approximation be valid, as dis- 
cussed above. There is also an upper bound on the range 
of spatial scales which comes about as follows. When the 
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instability is present, spatial Fourier modes with wave- 
length £ grow exponentially on a timescale 

^ (4.10) 




where is given by Eq. (|4.2p . If this timescale is longer 
than the Hubble time H^^, then the mode does not have 
time to grow and the instability is not relevant. Therefore 
the range of scales over which the instability operates is 



H 



(4.11) 



More generally, for a fluid of density p, if the instability 
is to be unmodified by the gravitational dynamics of the 
fluid, then the instability timescale must be shorter than 
the gravitational dynamical time, which from Eq. (|3.ip 



IS 



m-p / ^ PeS.{p) . Here Poff(p) is the total mass density 
of the coupled dark matter-dark energy fluid. This 
gives the criterion 



"T-off 



(4.12) 



which determines the values of density p and lengthscale 
C for which the instability operates. The upper length- 
scale can be rewritten using Eqs. (|4.6I1 . (|3.2|) and ()3.4p to 
give 



1 



™cft(p) 



™p|a'[0m(p)]| 

TOcff(/0) 



1 



1-- 



.13) 



Here the quantity d\aV/da(a) on the right hand side is 
expressed as a function of 4> using a = a{(j)), and then as 
a function of the density using </> = (pni (p) ■ 

At longer lengthscales, it is possible that the negative 
sound speed squared still engenders an instability, but 
determining this requires a stability analysis of the sys- 
tem in question including the effects of self gravity. The 
results might vary from one system to another. In this 
paper we will restrict attention to the regime (|4.12p where 
the presence of the instability can be easily diagnosed. 

The factor in square brackets in Eq. (|4.13p is always 
smaller than unity, since V is assumed to be a decreasing 
function of </> and hence also of a. It follows that the 
ratio of the maximum lengthscale £max to the minimum 
lengthscale £inin satisfies 

^<mp|a'[0„,(p)]|. (4.14) 

Hence in order for there to be a nonempty regime in 
which the instability operates, the strength of the cou- 
pling of the field to the dark matter must be much 
stronger than gravitational strength, 

mp|a'|>l, (4.15) 

as discussed in the introduction. In Sec. I IV CI below we 
give a simple explanation for this requirement. In Sec. 
IV Dl below we give an example of a model where we con- 
firm numerically instability is present at strong coupling 
but not when the coupling is weak. 



C. Jeans instability viewpoint 

There are two different ways of describing and un- 
derstanding the instability, depending on whether one 
thinks of the scalar-field mediated forces as being "grav- 
itational" forces or "pressure" forces. From one point of 
view, that of the Einstein frame description, the insta- 
bility is independent of gravity. This can be seen from 
the equation of motion (|3.ip : the instability is present 
even when the (Einstein-frame) metric perturbation due 
to the fluid can be neglected. In the adiabatic regime the 
acceleration due the scalar field is a gradient of a local 
function of the density [cf . Eq. (jAlOp below] , which can be 
thought of as a pressure. The net effect of the scalar in- 
teraction is to give a contribution to the specific enthalpy 
h(p) = J dp/ p of any fluid which is independent of the 
composition of the fluid. If net sound speed squared of 
the fluid is negative, then there exists an instability in 
accord with our usual hydrodynamic intuition. 

From another point of view, however, that of the Jor- 
dan frame description, the instability involves gravity. 
The gravitational force in this frame is mediated partly 
by a tensor interaction and partly by a scalar interaction. 
The effective Newton's constant describing the interac- 
tion of dark matter with itself is 



1 



2m2a'(</))2 
1 + ^ 



(4.16) 



where k is a spatial wavevector. This is Eq. (|Cip of 
Appendix [C] specialized to i = j = c, with Uc written 
just as a, and specialized to = (since experiments 
tell us that \a'^\ < 10~^m~^ today). Here the 1 in the 
square brackets describes the tensor interaction and the 
second term the scalar interaction. At long lengthscales, 
k <^ moff/(TOp|Q;'|), the scalar interaction is suppressed 
and we have Gcc ~ G. At short lengthscales, k ^ Woff, 
the scalar field is effectively massless and Gcc asymptotes 
to a constant, Gcc ~ G[l + 2mp(Q!')^]. However, when 
TOpla'l 3> 1 there is an intermediate range of lengthscales. 



mp\a'\ 



- < fc < rricft 



(4.17) 



in which the effective Newton's constant increases lin- 
early with A:^, 



Gc 



G- 



(4.18) 



'off 



A gravitational interaction with Gcc behaves just 

like a (negative) pressure in the hydrodynamic equations. 
This explains why the the effect of the scalar interaction 
can be thought of as either pressure or gravity in the 
range of scales (|4.17p . Note that the range of scales (|4.17p 
coincides with with the range (|4.13p derived above, up to 
a logarithmic correction factor. 

From this second, Jordan-frame point of view, the 
instability is simply a Jeans instability. In a cosmo- 
logical background with Einstein-frame metric ds^ = 
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—dt^+a{tY(lx^, the evolution equation for the CDM frac- 
tional density perturbation 5 with comoving wavenumber 
kr on subhorizon scales in the adiabatic limit is [33, ^3 



8 + 2H6 - 47rGcce"p(5 = 0, 



(4.19) 



where H = a/ a. Here Gcc is given by the expression 
(|4.16[) evaluated at the physical wavenumber k ~ kc/a, 
and we have neglected photons and baryons. 

Now in the absence of the Hubble damping term in Eq. 
(I4.19[) , the gravitational interaction described by the last 
term would cause an exponential growth of the mode, 
the usual Jeans instability of uniform fluid. Normally 
in a cosmological context, the Hubble damping term is 
present and the timescale ~ 1/H associated with this 
term is of the same order as the timescale 1 / \/Gp associ- 
ated with the gravitational interaction in the last term. 
Because of this equality of timcscales, the exponential 
growth is converted to power law growth by the Hub- 
ble damping. In the present context, however, things 
work differently. The gravitational constant governing 
the gravitational self-interaction of the mode is Gcc{k) 
instead of G, and consequently the timescale associated 
with the last term in Eq. (|4.19|) is shorter than the Hub- 
ble damping time by a factor of 



Gcc 

~G 



km-p\a'\ 



> 1. 



(4.20) 



Therefore the Hubble damping is ineffective and the 
Jeans instability causes approximate exponential growth 
rather than power law growth. 

The above discussion can also be cast in terms a scale- 
dependent sound speed instead of a scale dependent New- 
ton's constant; this clarifies the relation to our previous 
discussion of Sees. IIV Al and IIVBI As a slight general- 
ization, consider a fluid with an intrinsic sound speed 
Cs,in- Then the evolution equation (|4.19p generalizes to 
(see Appendix El 



where the effective total sound speed squared is 



c2. -l^Ge°p 



2mla'{(l)f 



k2 



(4.21) 



(4.22) 



In the range of lengthscales (|4.17p this squared sound 
speed is a constant, independent oik, as for a normal, hy- 
drodynamic sound speed. Outside of this range of scales, 
we have c^^t oc 1/fc at both large and small k (if the in- 
trinsic sound speed can be neglected), describing a con- 
ventional gravitational interaction. 

We reiterate that the existence of the range of scales 
(|4.17p in which Newton's constant scales linearly with fc^ 
depends on the assumption of strong coupling, |a'|mp <C 
1. If, instead, |a'|mp < 1, the dependence of G on k is 
very close to that of standard gravity, and the instability 
reduces to the normal Jeans instability of a fluid, charac- 
terized in a cosmological context by power law growth. 



D. Domain of validity of fluid description of dark 
matter 

Up till now we have described cold dark matter as a 
pressureless fluid. However, at a more fundamental level, 
one should use a kinetic theory description based on the 
coUisionless Boltzmann equation. In the conventional 
ACDM framework, the fluid approximation breaks down 
at small scales, below the free-streaming lengthscale, and 
also in the nonlinear regime after violent relaxation has 
taken place in CDM halos [1^. We now discuss how, in 
the models discussed here, the conventional picture for 
the fluid domain of validity is slightly modified. 

Let us denote by a the rms velocity of the dark mat- 
ter particles. Consider a perturbation characterized by 
a wavelength A and wavenumber k = 27r/A. The char- 
acteristic growth or oscillatio n time ass ociated with this 
perturbation is T(fc) ~ A/ \/|ctot(fc)p, where the total 
effective sound speed ctot is given by Eq. (I4.22p . The dis- 
tance traveled by a dark matter particle in this time is 
d{k) ^ aT{k), and the ratio of this distance to the size of 
the perturbation is 



d{k) aT{k) 



A 



(4.23) 



When this dimcnsionless ratio is of order unity or larger, 
perturbations do not have time to grow before they are 
washed out by free streaming of the particles, and the 
fluid approximation breaks down. Using the formula 
(|4.22p with Cs.in set to zero^^, , we obtain 



d{k) ak 
~ V47rGe"p 



2mla'{(t>f 



-I -1/2 



1 



k2 



(4.24) 



Now in the conventional CDM framework, the factor 
in the square brackets is unity, so the ratio (|4.24p is pro- 
portional to k and becomes large as fc — > oo. Hence the 
fluid approximation breaks down on small scales, below 
the critical free-streaming lengthscale Aps ^ oj^JGp. 

In the present context things work a little differently 
due to the scale dependence of Newton's constant. We 
can rewrite Eq. (j4.24|l in the approximate form 



i/2fcmp| OL 



d{k) 



/|c2| 



k < 



\/2mp|Q'| ' 



/|c2| moft 



nioff <C k 



(4.25) 

where = — 87rGe"pTOpa'^/m^jj is the constant value 
of the second term in the expression (|4.22p for c^^^. in the 



Since we expect Cs,in ~ o", setting Cg in to zero is only consistent 
in the regime <C |c^|. However, for > |c^|, retaining the 
intrinsic sound speed in Eq. I|4.23|l does not change the final 
result 114.261 1 for the free streaming scale in order of magnitude. 
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range of scales (|4.17|) . or equivalently the sound speed 
discussed in Sec. IIV Al We see that the ratio d[k)/\ is 
proportional to k at large scales and at small scales, but 
that in the intermediate range of scales it is a constant, so 
that the effect of free streaming is equally important for 
all the modes in this range. If we define the free stream- 
ing Icngthscalc AfsI^) to be the smallest Icngthscale for 
which free streaming is unimportant, d{k)/\ < 1, then 
we obtain 



Afs(ct) 



1 



/|c2| '"ott 



1^ < 



(4.26) 



This lengthscale jumps discontinuously at ^ jc^l. 
There are thus two different regimes that occur: 

• When (7^ ^ |c^|, free streaming is important only 
at scales small compared to l/mes (for which the 
adiabatic approximation is invalid anyway). The 
fluid approximation is valid throughout the range 
of lengthscales (|4.17p . and so the adiabatic insta- 
bility is present. This conclusion is confirmed by 
a kinetic theory analysis (see Appendix [Bjl. which 
shows that linearized perturbations of any homo- 
geneous, isotropic initial particle distribution func- 
tion are always unstable on scales that are in the 
regime (|4.17p . as long as tr^ ^ \cl\. 

• When > |cg|, free streaming becomes important 
and the fluid approximation breaks down through- 
out the range of scales (|4.17p . One expects the free 
streaming (also called Landau damping) to kill the 
instability. This is confirmed by our kinetic the- 
ory analysis of Appendix [B| we show that for a 
Maxwellian distribution, the finite velocity disper- 
sion stabilizes the coupled fluid whenever cr^ > |c 
in agreement with the analysis of Ref . [lB| • 
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Consider now the evolution of cosmological perturba- 
tion modes. When will our analysis of instability apply? 
First, the condition < |c^| is not very restrictive, since 
the CDM cools rapidly with the Universe's expansion, 
with temperature scaling as (1 + z)^ once the particles 
become non-relativistic. For example, for CDM particles 
of mass ~ GeV with weak scale cross sections, the CDM 
temperature is ^ 10"'' K at decoupling [s^. However, 
after perturbations go nonlinear and violent relaxation 
takes place in CDM halos, the effective coarse-grained ve- 
locity dispersion becomes much larger. Hence, our anal- 
ysis does not apply to modes that arc in the nonlinear 
regime. Our analysis will apply in the early Universe, 
before any modes have gone nonlinear. It will also ap- 
ply to large scale modes, even after smaller scale modes 
have gone nonlinear, since such large scale modes should 
still be well described by linear theory (see for example 
the qualitative arguments in chapter 28 of Peebles [3q). 
Our investigation of specific models later in the paper 
will focus on these large scale, linear regime modes. 



We note that earlier investigations of the instability 
focused instead on small scale modes, below the free 
streaming scale Aps [H, [III • From the formula (|4.26p we 
see that, for the models discussed here, either the adi- 
abatic approximation is not valid on these small scales 
since Aps ^ ^/nicS, or > \cj.\ and the instability is 
killed by free streaming. 



V. EXAMPLES OF THEORIES WITH 
ADIABATIC INSTABILITY 

In this section we discuss some specific classes of the- 
ories. 



A. Exponential potential and constant coupling 

We first consider theories with exponential potentials 
of the form 

V = l/oe-^'^/™p , (5.1) 
with A > and with linear coupling functions 

a{(j)) = -pC— , (5.2) 

TOp 

where 13 = •\/2/3 and C is a constant. These theories 
have been previously studied in Ref. Q. The effective 
potential is, from Eq. (|2.1ip . 



(5.3) 



and solving for the local minimum of this potential yields 
the relation between (f> and p in the adiabatic regime: 



,(A-/3C)0„,(p)/mp _ 



-l3Cp 



(5.4) 



Note that C must be negative in order for the effective 
potential to have a local minimum and for an adiabatic 
regime to exist. We will restrict attention to this case, 
and we define the dimensionless positive parameter 7 = 
— A//3C. The corresponding effective mass parameter is 



'off 



A'm-Vo 



1+7 f P 
7 VTK) 



(5.5) 



Next, we compute the sound speed squared. Using Eq. 
(14.21) we obtain 



1 



1 + 7 



(5.6) 



^'^ The notation in Eq. I|5.2| l is chosen such that f{K) gravity theo- 
ries correspond to C = 1/2. 
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so this model is always unstable in the adiabatic regime. 
From Eqs. (|5.ip and (|5.4p we also obtain 



7 



ainy 



din p 1+7 



(5.7) 



We now insert the effective mass (|5.5|) . the sound speed 
squared (|5.6p and the logarithmic derivative (|5.7|) into 
Eq. (|37T5)) and into the second half of Eq. ((47T2)) . This 
yields the range of spatial scales £ over which the insta- 
bility operates for a given density p to be Cmin{p) ^ £ <C 
CrmiAp)i where 



7 



A2(l+7)' ^0 V P 



and >Cmax(p)^ is a constant times this: 



(5.8) 



(5.9) 



Thus, there is a nonempty unstable regime only when 
/3|C| ^ 1, ie with the scalar coupling is strong compared 
to the gravitational coupling, in agreement with the dis- 
cussion in Sec. IIVBI 

To see the effect of the instability more explicitly, we 
consider cosmological perturbations. The Einstein- frame 
FRW equation in the adiabatic limit is 



V ^ 



(5.10) 



where p oc 1/a^. This yields a{t) cx i2/(3+3u)„,t)^ where 
the effective equation of state parameter is 



1 



(5.11) 



In the strong coupling limit |C| ^ oo that we specialize 
to here, —1. Thus the adiabatic regime of this 

model with large \C\ is incompatible with observations 
in the matter dominated era, where we know Woff ~ 
except for at small redshifts. Nevertheless, the model is 
still useful as an illustration of the instability. 

We find from Eqs. (jO)) . ([g^ and ((ETU)) that the range 
of unstable scales is given by 



1 



/32C2 



1 



fc2 3(1 + 7) 



(5.12) 



where k is comoving wavenumber. This range of scales 
always lies just inside the horizon. A given mode k will 
evolve through this unstable region before it exits the 
horizon. 

Next we use the approximate form (j4.18p of Newton's 
constant in the perturbation evolution equation (|4.19p . 
and transform from t derivatives to a derivatives. This 
gives 



cP^ ^ 3 / I rilnpcff 
da? a \ 2 In p 



d5 
da 



(gQ^fc^e^p 
ml^H'^a-^ 



(5 = 0. 
(5.13) 



Specializing this equation to the exponential model using 
Eqs. dEH), ([53)1 and ((STTU)) and taking the strong 

coupling limit |C| oo gives 



da"^ 



3 dS 
a da 



H'^a^ 



(5.14) 



In the strong coupling limit H is approximately a con- 
stant, H « Hq, and the growing mode solution is 



a \-tiQa 



ttHo 
2ka 



exp 



k 



Hoa 



(5.15) 



where Ki is the modified Bessel function. The mode 
grows by a factor ~ e when the scale factor changes from 
a to a -|- Aa, where Aa/a ^ aHo/k <^ 1 for subhorizon 
modes. 

A more detailed analysis of the cosmology of this model 
is given in Ref. [Ill, but in the non-adiabatic regime 
|C| ~ 1 rather than the strong coupling regime |C| » 1 
considered here. 



B. Two component dark matter models 

We next consider models in which there are two dark 
matter sectors, a density pc which is not coupled to the 
scalar field, and a density pco which is coupled with cou- 
pling function (|5.2p and exponential potential (|5.1|) . Both 
of these components are treated as pressureless fluids. 
The FRW equation for this model in the adiabatic limit 
is [cf. Eq. (|5l0l) above] 



3mlH^ + e^pco + Pc 



(5.16) 



Similar two component models have been considered by 
Farrar and Peebles [2^ . This model is also similar to the 
mass varying neutrino model model p^ . [iTj where the 
neutrinos play the role of the coupled component; see 
Sec. IV Dl below. The first two terms on the right hand 
side of Eq. (|5.16[) act like a fluid with equation of state 
parameter given by (|5.1ip . and in the strong coupling 
limit |C| ^ 1 this fluid acts like a cosmological constant. 
Thus, the background cosmology can be made close to 
ACDM by taking \C\ to be large. 

The fraction of dark matter which is coupled must be 
small in the limit of large coupling, |C| ^ 1. Denoting 
flv = V/{3mlH^), i},o = e"p™/(3m2ir2) and fl, = 
Pc/ (SmpH^), we have 1 = ily + ^co + ^c- Also from Eq. 
(|5.4p it follows that, if the asymptotic adiabatic regime 
has been reached, flco = J^v, and we we obtain 



H-7 



(5.17) 



Since ilc ^ 0.3 today, and 7 <C 1 in the strong coupling 
limit we are considering, we must have f2co *C 1 today. 

The maximum and minimum lengthscales for the in- 
stability are still given by Eqs. (|5.8p and (|5.9p . but with 
p replaced by pco- Since pco is approximately a constant 
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in the strong coupling limit, these lengthscales are also 
constants. If the parameters of the model are chosen so 
that ^Ic ^ ^ today, then 



(5.18) 



The evolution equations for the fractional density pertur- 
bations (5j = ^Pj/Pj in the adiabatic limit on subhorizon 
scales are given by 



(5j + 2HSj - 47r^ GjkPke"'=(5k = 



(5.19) 



where the effective Newton's constants are given by Eq. 
(|C1[) with ab set to zero. Writing this out explicitly we 
obtain 



(5.20) 



(5co + 2H5co = 



1 



2ml 



PcSc 



2ml 



e"PcoSc 



(5.21) 



The condition for the instability to operate is that the 
timescale associated with the second term on the right 
hand side of Eq. (|5.2ip be short compared with H~^, or 



(5.22) 



Now the effective mass for this model is given by m^^ = 
(3^C^m-'^{V + e"pco) - 3(3^C^m~^H^{nv + i^co)- Sub- 
stituting this into Eq. (|5.22[) and using Eq. (|5.17|) gives 
the criterion k/{aH) ^ 1. Therefore the instability 
should operate whenever modes are inside the horizon 
and in the range of scales (|5.18p . 

These expectations are confirmed by numerical inte- 
grations. In figure [1] we present a numerical analysis 
of such a two component model. We consider an ex- 
ponential potential with A = 2 and strong coupling with 
C = —20, and typical cosmological parameters are as- 
sumed. Ho = 70kms-iMpc"\ n,, = 0.05, Qc = 0.2, 
fico = 0.05, and fiy = 0.7. We fix initial condi- 
tions of 4)/m-p = 10^'^ and ^ = at a = 10^^°, al- 
though the dynamical attractor renders the final evolu- 
tion largely insensitive to these choices, we have checked 
that (j)/mp{a = 10~^°) = 10~^° — 1 give the same evo- 
lution. We neglect the effect of the coupled CDM com- 
ponent peculiar velocity in the initial conditions (it is 
many orders of magnitude smaller than the density per- 
turbation) and assume that the coupled and uncoupled 
CDM components have the same initial fractional den- 
sity perturbations 6c = 5co, fixed by the usual adia- 
batic initial conditions. As shown in figure [U the back- 
ground evolution is entirely consistent with a ACDM 
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FIG. 1: [Bottom] The two component coupled dark energy 
(CDE) model, with exponential potential and coupling, with 
A = 2 and coupling C = -20 with Ho = 70kms~^ Mpc"\ 
Q.b = 0.05, Qc = 0.2, Qco = 0.05, and Qv = 0.70. At 
late times the scalar field finds the adiabatic minimum with 
asymptotic equation of state, and sound speed = — l/(l-f7) = 
—0.89, able to reproduce a viable background evolution con- 
sistent with supernovae, CMB angular diameter distance and 
BEN expansion history constraints. The figure shows the evo- 
lution of the eff'ective equation of state, Weff — Ptot/ptot = 
(2/3)(dlnf/dlna) — 1, (black full line), the adiabatic speed 



of sound. 



P/p for all components (blue long dashed 



line) and for the coupled components only (green dot long 
dashed line), and effective speed of sound for = 5P/5p 
at = 0.01/Mpc for all components (red dot-dashed line) 
and for the coupled components alone (magenta dotted line). 
The efi'ective equation of state for a comparable ACDM model 
with ric = 0.25, Q,b = 0.05 and Q.a = 0.7 is also shown (black 
dashed line). [Top] The growth of the fractional over-density 
5 — 5p/p for k = 0.01/Mpc for the coupled CDM compo- 
nent, Sco, (red long dashed line) and uncoupled component, 
Sc, (black full line) in comparison to the growth for the ACDM 
model (black dashed line). At late times the adiabatic behav- 
ior triggers a dramatic increase in the rate of growth of both 
uncoupled and coupled components, leading to structure pre- 
dictions inconsistent with observations. 



like scenario, with Wcff = —0.69 today and an asymp- 
totic equation of state at late times, given by (|5.1ip . 
Wcff = —1/(1-1-7) = —0.89. The large coupling drives the 
evolution to an adiabatic regime at late times, with an 
adiabatic sound speed — s- —1/(1 + 7) as in (|5.6p . This 
drives a rapid growth in over-densities once in the adia- 
batic regime so that although consistent with structure 
observations at early times, they arc inconsistent once 
the accelerative regime has begun. 
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In summary, these models provide a class of theories 
for which the background cosmology is compatible with 
observations, but which are ruled out by the adiabatic 
instability of the perturbations. 



C. Chameleon models 

Next we study the so-called chameleon models [l^ [2C 
defined by the potential 



= AM* ( ^ 



(5.23) 



where M is a mass scale and n > and A are dimension- 
less constants, together with the coupling function (|5.2[1 . 
In these models it has been previously shown that the 
adiabatic regime is achieved in static solutions describing 
macroscopic bodies like the Earth, and that cosmologi- 
cal solutions in the adiabatic regime provide good models 
of dark energy [13, SOj EH ■ We now study under what 
conditions these models are unstable. 

The effective potential is, from Eq. (|2.1ip . 



(5.24) 



and solving for the local minimum of this potential yields 
the relation between <f> and p in the adiabatic regime: 



P 



(5.25) 



Here x — —(3C4>c{p)/mp is dimensionless and the critical 
density is 



Pcrit = nXM 



nip J 



(5.26) 



As before the existence of a local minimum in the effective 
potential requires C to be negative. We shall restrict 
attention to the regime 



P > Pcrit 



(5.27) 



since for models of dark energy pcrit will be of order the 
present day cosmological density. In this regime the so- 
lution to Eq. (|5.25p is approximately 



The corresponding effective mass parameter is 
.2 .^.^,.2^~/^CM^"+^ 



is = {x + n + l)nXAr 



xnip ) 



which in the regime (|5.27p simplifies to 



= (n + l)(-/3C)2mpVrit 



Pcrit 



2+2 

T+1 



(5.28) 



(5.29) 



(5.30) 



Next, we compute the sound speed squared. Using Eq. 
(|4.2p we obtain 



1 



-1 



1 TOc 



(-/3C) <t> ' 



(5.31) 



and since C is negative, we see that this model is always 
unstable in the adiabatic regime. Inserting the effective 
mass (|5.29p and the sound speed squared (|5.3ip into Eq. 
(|4.12p we obtain the range of spatial scales £ over which 
the instability operates for a given density p: 



CrniApf < « {PCfl^^^AP? 



where 

^mi„(p) = 



Pcrit 

(7l+l)(/3C)2p„,t V~ 



(5.32) 



(5.33) 



We see that the range of unstable lengthscales is non- 
empty only if 



» 1, 



(5.34) 



which as before is equivalent to the strong coupling con- 
dition (|4.15p . Note that the first of the two inequalities 
in Eq. (15.321) is equivalent to the "thin-shell condition" 
of Ref. [201 when the background value of the scalar field 
can be neglected. 

As for the exponential models of Sec. IV A[ the effec- 
tive equation of state is close to — 1 in the adiabatic 
regime for large coupling, assuming p ^ Pcrit- Therefore 
these models do not give an acceptable background cos- 
mology for the matter dominated era in their adiabatic 
regime. However, one can construct two component mod- 
els analogous to those in Sec. IV Bl using the chameleon 
potential. Those models give an acceptable background 
cosmology, but are then ruled out by the adiabatic insta- 
bility. 



D. Mass Varying Neutrino ("MaVaN") models 

The impact of adiabatic instabilities has been discussed 
extensively in the context of MaVaN models, in which 
the light mass of the neutrino and the recent accelerative 
era are twinned together through a scalar field coupling 
[12, [Hj [13 • The adiabatic instability was shown to be a 
concern in these models with the implication of forming 
compact localized regions of neutrinos after undergoing 
dramatic adiabatic collapse [Tsj. 

The action for the MaVaN models is of the form (|2.2p 
with p = 0. Therefore the analysis in the earlier sections 
of this paper, which assumed p = 2, does not directly 
apply to these models. However it is straightforward to 
generalize our analysis to cover this case. Consistent with 
the general discussion in sections IIIHIIVI the instability 
is present in these theories if the coupling is strong. Here 
we discuss examples of models which evade or are subject 
to the instability on cosmological scales. 
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Recently, Ref. [17| discussed a MaVaN scenario with a 
logarithmic potential and scalar dependent mass, 



Folog(l+e0) 



(5.35) 
(5.36) 



where ^ is a constant and m^o and (f>^ are the current 
values of the neutrino mass and the scalar field respec- 
tively. It was found in p7j that this model can exhibit 
an instability in growth for an otherwise cosmologically 
viable background solution. We find, however, that this 
model can also allow stable solutions for identical frac- 
tional densities today as those studied in [l3| , for a wide 
range of parameter values. If the coupling, m'^/m^, is 
not large compared to m^^ at late times, the evolution 
never enters the adiabatic regime on cosmic scales. This 
translates in this model to (j)s, not being significantly less 
than TOp. 

The MaVaN model evolves according to a coupled 
Klein Gordon equation. 

4> + 2H^ + a^V'{^) = -a^ipU^)-3P;,{^)) . (5.37) 

Assuming no chemical potential, and g spin states per 
neutrino species with momentum p and mass m ,the neu- 
trino density and pressure are given by 



a Pu = 



27r2 



dqq\q^+a^ml)^f{q) (5.38) 







„4p _ giksTsr r,„„2 

Qj 1,, - 



27r2 

fiq) « [exp(rj) + l]-^ 



dqq 

3(g2 + a^TO^)? 



r/(g)(5.39) 
(5.40) 



with q = ap/ksT^ and m,y = m^c^ /ksT^ ■ 

In the relativistic regime, with m ^ 1, the potential is 
negligible and the driving term, on the right hand side of 
(|5.37p , can be calculated by doing a Taylor expansion to 
first order in m^, 

a\p. - 3P.) « I dqqf{q)a'ml (5.41) 



so that 



a'{pl-3Pi) 



-ip.-3P,)a^ 



(5.42) 



(5.43) 



where po = Ttt 3(^5 T") / 240 would be the relativistic 
neutrino energy density per neutrino species today with 
temperature T^. 

Putting (|5.43p into (|5.37p and neglecting the potential 
we find a power law attractor </> oc with x = 0.5. The 
normalization of 4> is wholly specified in the attractor 
by (inSZl)- Writing </> = 0j(r/Ti)"-^, and a oc tP , with 
p = 2/(1 + 3wcff), we find 



SOpoml^cfl 



77r2(4|j- 1) 



^0.5 



(5.44) 
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FIG. 2: [Top panels] Evolution of m^((/)) and neutrino tem- 
perature (left), and associated equation of state and adiabatic 
sound speed (right) for MaVaN model described in the text 
with m^o ~ 0.312el/ and (^t = 1.8mp. [Bottom left panel] 
Scalar field evolution in the MaVaN scenario, for 3 values of 
(j>t ~ 10~^mp, 0.3mp and l.Smp (full lines) showing the r"'^ 
attractor while the neutrino is relativistic allowing late time 
evolution to be independent of initial conditions. The vac- 
uum expectation value of the scalar field, if the field becomes 
adiabatic, (dashed lines) is also shown. For > 10~^mp the 
scalar field does not enter an adiabatic era on cosmological 
scales before now, and growth of perturbations remains well- 
behaved. For smaller (^*, for example 0, ~ 10~^mp shown, 
the evolution is adiabatic at late times, similar to that dis- 
cussed in [13] • [Bottom right panel] The resulting matter 
power spectrum from the coupled dark energy (CDE) model 
with (j>, = l.Smp and = 0.02 (full line) is very similar 
to that for ACDM (dashed line) with the same baryon frac- 
tion and Ho, fi^ = 0.3, fii^ = 0.02 when normalized at large 
scales. 



When the neutrino is non-relativistic, if we again ne- 
glect the potential, 



(p, -3P,) = 
a'{pl^3P;,) = 



(5.45) 
(5.46) 



The Klein-Gordon equation has a solution (j) with 
x = (2 — p)/3, tending towards a cessation of growth in 
<j) in the matter dominated era. 

If the neutrinos are non-relativistic (but, of course, as- 
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suming that they decoupled when they were relativistic) 

(5.47) 



180C(3) po 
77r4 a3 



1/0 



(5.48) 



At late times, the complete effective potential is relevant, 
for which there exists a minimum at positive (j) given by 



1 VL 



a3 n 



^iog(i + e0o . 



(5.49) 



pot 



We modified CAMB [39[ to investigate the evolution 
numerically. In figure [2] we show the scalar field and 
neutrino mass evolution for ^ = lO^'^mp^, (j)^, = l.Smp, 
and m^o = 0.312ey (giving = 0.02), ^i, = 0.05, 
= 0.23, Hq = 70kms~iMpc"^ In the bottom left 
hand figure, the numerical evolution of scalar field is 
shown for 0* = 1.8mp along with two smaller values 
0* ~ 10~'^TOp and O.STOp for which the minimum of the 
effective potential is steeper. The background evolution 
obeys the ^ r°-^ attractor and normalization in (|5.44[1 . 
rendering it largely independent of the initial conditions. 
When the neutrino becomes non-relativistic the evolution 
slows and for small values of (/)* starts to track the VEV, 
as discussed in [l3]. However for larger values of <^h., the 
VEV is not reached until later times, and in addition, the 
effective potential minimum is shallow enough that the 
scalar field does not get fixed at the VEV automatically, 
enabling well-behaved growth. The bottom right hand 
figure shows the resultant matter power spectrum, which 
is very similar to a fiducial ACDM model with the same 
Hubble factor, baryon and neutrino density today when 
normalized to the same amplitude at large scales. 



VI. CONCLUSIONS 

If dark energy and dark matter are to fit into a coher- 
ent fundamental physics framework then there are likely 
couplings between them. Such couplings may have far 
reaching macroscopic implications on scales ranging from 
the solar system up to cosmological horizon. As we have 
discussed, these effects may lead to observationally dis- 
tinctive characteristics, which may allow us to tease out 
the nature of the dark sector. However they may also 
give rise to catastrophic instabilities with which we may 
constrain the class of physically viable dark energy mod- 
els. 

In this paper wc have considered such theories in which 
there exists a nontrivial coupling between the dark mat- 
ter sector and the sector responsible for the acceleration 
of the universe. 

We have comprehensively analyzed an instability - 
characterized by a negative sound speed squared of an 
effective coupled dark matter/dark energy fluid - that 
exists whenever such theories enter an adiabatic regime 



in which the scalar field faithfully tracks the minimum 
of the effective potential, and the coupling strength is 
strong compared to gravitational strength. The adia- 
batic regime occurs when the relaxation time scale as- 
sociated with the scalar field is much shorter than the 
Hubble time. We have demonstrated how this instability 
can be viewed from the alternative perspectives of the 
kinetic theory of dark matter and as a Jeans instability 
associated with modified Newton's constants. 

We have established the conditions under which the 
adiabatic instability occurs, finding a condition on the 



coupling. 



^ 1, and have identified the time and 



length scales over which the instability is active in a given 
setting governed by the matter density. These length 
scales can differ greatly dependent on whether one is con- 
sidering galactic or cosmic densities. 

Our work builds on previous analyses of MaVaN 
[H, [13, [13, m, El , chameleon [ii,!!^, and general cou- 
pled scenarios [l6|. Our numerical analyses of coupled 
CDM and MaVaN models reinforce our analytic find- 
ings. We show that the regime of adiabatic behavior 
agrees with the predictions based on coupling strength 
and length scale conditions mentioned above. In the ap- 
propriate limits, our results reduce to the previous find- 
ings in a number of cases, while in several other cases we 
have provided corrected results. In particular, we have 
shown that stable MaVaN models exist that evade the 
instability, if m'^jmi, is not large compared to Trip. 
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APPENDIX A: EFFECT OF NORMAL MATTER 
ON INSTABILITY 

In the body of this paper we have neglected the effect of 
baryons and the other matter species, assuming that their 
density is much smaller than that of the dark matter. A 
natural question is does the instability persist in regimes 
where the density of visible matter density is comparable 
to or larger than that of the dark matter - one might 
expect the instability to be killed by the pressure of the 
visible matter. In general the answer depends on the 
coupling function a\y in the action (|2.ip . which governs 
the coupling of the scalar field (/) to normal matter. The 
instability persists in two specific cases, (i) a;, = where 
normal matter is uncoupled from and (ii) a\, = etc, the 
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case usually considered in which dark matter and normal 
matter couple the same way to the scalar field, respecting 
the equivalence principle. 

To derive this result, we generalize the derivation given 
in the body of the paper from a single, pressureless fluid 
to a set of N non-interacting fluids, each of which can 
have a pressure. We start from the action (|2.ip . and we 
define the Jordan frame metric that the jth sector couples 
minimally to by 



gjab 



c2aj(0) 



9ab 



(Al) 



The stress energy tensor T-"-'' of the jth sector is defined 
by 



Sj [g, ab+Sgj ab, '^j] = Sj [gj ab, + ^ 



d^3 



^g;T.;'^5-g;ab. 
(A2) 

and we assume it has the form of a perfect fluid 



ah 



ah 



(A3) 



Here is the 4- velocity which is normalized according to 
gjabUj^uf' = —1. This perfect fluid assumption requires 
that the matter in the jth sector be barotropic, ie that 
its pressure be determined uniquely by its density, which 
in many regimes is a good approximation. 
The equation of motion of the jth fluid is 



(A4) 



where Vj a is the derivative operator determined by the 
metric gjab- The components of this equation perpendic- 
ular to and parallel to the 4-velocity are 



1 



Pi +Pj 



-ig^'' + ufu^)VibPi 



(A5) 



where af = uf'W\ huf is the Jordan frame 4-acceleration. 
and 



(A6) 



Next, we conformally transform these fluid equations us- 
ing Eq. (|A1[) to write them in terms of the Einstein frame 
metric gab and the Einstein-frame normalized 4- velocities 



e^'u^ which satisfy gabu^u^ = — 1. This gives 



Pi +Pi 



where 



u^Wbuf^ and 

-3a; w r^3a 



e-3"^V, [e^"'{pj = ufVaP^ 



(A7) 



(A8) 



Next, we consider the non-relativistic limit of the fluid 
equations (jA7[) and (|A8|) together with the adiabatic limit 
of the scalar field equation (|2.4p . We also neglect self- 
gravity in the Einstein frame, taking gab ~ Vab- This 



approximation should be valid on sufficiently small spa- 
tial scales. We get 



d_ 
di 







(9v 
'dt 



^+(vj.V)vj=-^-Vaj, 



Pi 



(A9) 



(AlO) 



where aj = and = (t>{pk) is given by the equation 

r(0) = -^a;(</,)e4".Wp. . (All) 
j 

Next, we switch to using the rescaled density variables 

Pi = e'"^Pj . (A12) 

Wc linearize the resulting equations about the back- 
ground solution Pj = Pjo = constant, vj = 0, = 0o- 
In other words we assume that the dark matter and vis- 
ible matter densities are constant and that the fluids 
are not in relative motion. We define ajo = Q;j((/)o) and 
= a'-^{(j)o), and look for a solution 



Pi = PiO + ^Pi^ 



(A13) 
(AM) 



From Eq. (jAll[) we obtain that 54) = Xj<^Pji where 

(A15) 



Xj 



V" 



EkKo + Ko)'] 



'PkO 



Next we assume that all the variables are proportional 
to exp [ik • X — iujt] , and that the velocities are propor- 
tional to k. This assumption will yield one particular set 
of modes of oscillation of the coupled fluids, but these 
modes will contain the instability if it is present. The 
resulting eigenvalue equation for uj'^ is 



(A16) 



where the matrix Fji;, which plays the role of the effective 
squared sound speed, is 



Here 



Tjk = Cj^(5jk + (1 - 3cj^)/3joajoXk 



dPi 



dPi 



(A17) 



(A18) 



is the physical (Jordan frame) squared sound speed of 
the jth fluid. 

From the eigenvalue equation (|A16p , it follows that the 
system will be stable if and only if all of the eigenvalues 
of the matrix Fjk are real and positive. For the case of 
a single fluid with c? = 0, this criterion reduces to the 



criterion (|4.3p derived in the body of the paper. In gen- 
eral, there is a competition between the positive squared 
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sound speeds of the fluids in the first term in Eq. (jA17p . 
and the (possibly) negative squared sound speeds com- 
ing from the interaction with the scalar field given in the 
second term. 

We now specialize to two fluids, a dark matter fluid 
with zero sound speed, and a fluid describing the visible 
matter. Taking j = 1 = c for the dark matter (CDM) and 
j = 2 = 6 for the visible matter (baryons), and defining 

v; = {l- 3c2)pjoa;o (A19) 
gives for the matrix 



jk 



VcXc 
VbXc 



^cXb 



VbXb 



(A20) 



Now if the system is stable, then both of the eigenval- 
ues Ai and A2 of the matrix must be real and nonnega- 
tivc, and hence both the determinant A1A2 and the trace 
Ai -|- A2 of the matrix must be nonnegative. Conversely, 
if either the determinant or the trace of the matrix is 
negative, then the system is unstable. The determinant 
is 



detr 



CbVcXc 



K) 



4pc 



(A21) 



where 



+ [a'l + [a'.f] e"V6 (A22) 

and where we have used Eqs. (jAlSp and (|A19p . In Eqs. 
(IA21[) and (jA22p , can be taken to be the function of pc 
and given by [cf. Eq. (|A8p ] 

V'{(f) = -a',(0)e"=(*Vc - aU0)e"'^*V^ • (A23) 

Equations (jA2ip - (|A23p allow us to determine which val- 
ues of Pc and pb satisfy the sufficient condition det F < 
for instability, given the functions V{(l)), ac(4>) and Q;f,(0). 
Note that ac('/>) was denoted by a{(t)) in the body of the 
paper. 

We now consider two special cases. If = 0, so that 
the visible matter is not coupled to the scalar field, then 
eigenvalues of the matrix T arc just and VcXc- The 
second of these eigenvalues coincides with the effective 
squared sound speed computed above in Eq. (|4.2[) . Thus 
we recover the results of the body of the paper. 

The second special case is when = ac- In this case 
the instability criterion det F < again reduces to the 
criterion (|4.2[) computed earlier for a single fluid, with 
the modification that the effective sound speed squared is 
now a function of the total density p = Pc+Pb rather than 
just of Pc- If this total density is in the unstable regime, 
then the instability will persist despite the pressure of 
the visible matter. 

Finally, one can also consider the effect of the scalar 
field on just the baryons, neglecting the dark matter. 
The sound speed of the baryons gets an additional term, 
but it is a small correction unless the coupling is large, 
aJ^TOp ^ 1, and observational tests of general relativity 
in the Solar System require aj^mp <^ 1. 



APPENDIX B: KINETIC THEORY TREATMENT 
OF INSTABILITY 

In this appendix we describe dark matter using the col- 
lisionless Boltzmann equation, and specialize to the non- 
relativistic regime and to lengthscales small enough that 
self-gravity can be neglected. We show, first, that the 
instability is generic, occurring for any velocity distribu- 
tion function with sufficiently small velocity dispersion, 
and; second, that for MaxwcUian distributions the fluid is 
stabilized by free streaming once the velocity dispersion 
becomes larger than a critical value. 

Starting from the general action (|2.ip . we special- 
ize to the non-relativistic limit neglecting self-gravity 
{gab ~ Vab), which wiU be valid on sufficiently small spa- 
tial scales. The Jordan-frame metric is then 

We assume that dark matter is composed of non- 
interacting, non-relativistic particles of mass p. We de- 
note by /(t,x, v) the one-particle distribution function, 
normalized so that the number of particles in the volume 
element d^x and in the velocity region d^v is 

f(t,x,\r)d^xd^v . 

Now the physical (Jordan-frame) volume element is 
e^^d^x, so the Jordan frame mass density p [cf. Eq. (|2.3[) 
above] is given by p = pe~^°' J d^vf. From Eq. (|2.7p the 
rescalcd density variable p is then given by 



p^ p / d vf 



(Bl) 



With this notation the coUisionless Boltzmann equation 
is [cf. Eq. (jAlOp above with pj = 0] 



dl 
dt 



.9/ 
dx' 



da df 

dx^ dv^ 



(B2) 



In the adiabatic limit, the variable a in this equation is 
given by a = «(</)), where (j) is given in terms of p by Eq. 
(|3.4p . and p is given in turn in terms of / by Eq. (jBip . 

We now linearize the Boltzmann equation about a ho- 
mogeneous background solution, taking / = /o(v) + 
(5/(t,x, v). and taking 8f oc exp[— zwt -I- ik • x]. Using 
the identity da/ dp = c^/p from Eq. (|4.2p . and integrat- 
ing over V, gives the dispersion relation 



(B3) 



Here po is the background density. Without loss of gen- 
erality we take k to be in the z-direction, assuming /o is 
isotropic, and we define 



dVr, 



dVy fo{Vx,Vy,v) 



(B4) 
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This distribution function is normaUzed so that J dvF 
1, and the dispersion relation can now be written as 



1 + c 



dv 



F'{v) 
Lo/k — 



(B5) 



Given the distribution function F{v) we can solve this 
equation to obtain to = ti'(fc), which is complex in general. 

The integrand in Eq. (jBSP contains a singularity at 
V = uj/k. The derivation given here is incomplete since 
it does not provide a specification for how to deal with 
the singularity. However, just as in plasma physics [i^ . 
it is possible to give an alternative derivation based on 
Laplace transforms. That alternative derivation yields 
the following specification: the integral over v must be 
taken over the Landau contour, which runs along the real 
axis if Im(aj) > 0, but dips below the real axis to encircle 
the pole at w = (jj/fc if Im(w) < 0. 

We now specialize to the adiabatic regime where < 
0, and we write — ~0^. We write the complex fre- 
quency in terms of its real and imaginary parts, uo = 
LJr + ioji, and we look for a solution w = uj{k) of the dis- 
persion relation (jB5[) with uji > 0, corresponding to an 
unstable mode. For such a solution the Landau contour 
is along the real axis which simplifies the analysis. 

If the velocity dispersion of the distribution F{v) is 
small, we can expand the denominator of the integrand 
in Eq. (|B5[) as a power series in vk/cu. Integrating by 
parts and solving the resulting equation for cj^ gives 



UJ 



-0^ + 3a^ + 



(B6) 



where cr^ = / dvv'^F{v). This equation has a solution 
with positive imaginary part, > 0, consistent with 
the assumption used in its derivation. Therefore the in- 
stability is generic, present for any velocity distribution, 
as long as CT ^ /9. Equation (|B6p also suggests that 
the instability will be removed when a gets large, as ar- 
gued in Ref. p^ . since then the positive second term in 
the square brackets will overcome the negative first term. 
However, Eq. (jB6p is only valid in the regime <t <C /3, and 
so to investigate stabilization we must use an alternative 
method of computation. 

Returning to the general dispersion relation (jBSp . we 
obtain from its real and imaginary parts the equations 



and 



d., .=0 



, (B7) 



(B8) 



[oJr ~ kvY + uf 

We now specialize to the Maxwellian velocity distribution 



F{v) 



2TTa 



-e 2<T 



(B9) 



For this case Eq. (|B8|) has a unique solution for tUr, 
namely tOr — 0. This can be seen by splitting the inte- 
gral into V > and v < contributions and substituting 
V — > —V in the w < term. This yields 



= LOr 



dv 



[{Ujr - kvY + Ujf] [{uJr + kvf + Uj'j 



(BIO) 

and the result follows since the integrand is everywhere 
positive. Equation (|B7|) now simplifies to 



dv 



V^ -\- 



1 



/27r 



(BU) 



where we have defined v = v /a and n = LOi/ (ka). 

Equation (|B11[) determines k and thence oji as a func- 
tion of cr//?. When k is large, expanding the factor in 
round brackets in the integrand as a power series in v/k 
gives the result (|B6|) above. As k decreases, a/P in- 
creases, until as K 0, a/ (3 approaches the limiting value 
a//3 = 1. It can be seen that there are no solutions to 
Eq. (jBlip with a/P > 1, since the right hand side is a 
monotonic function of k^. Therefore, whenever^^ 



a> P 



(B12) 



there are no unstable modes with uJi > 0. In other words, 
the instability has been removed by the damping process 
associated with the finite dispersion (Landau damping, 
also called free streaming). The perturbation to the dis- 
tribution function is proportional to the integrand in Eq. 
(|B11[) . Therefore the nature of the damping is that in- 
teraction with the growing mode moves some particles 
from velocities w ~ /? to larger velocities, removing en- 
ergy from the mode. 



APPENDIX C: EFFECTIVE NEWTON'S 
CONSTANT 

In this appendix we derive the formula (jCip for the 
effective Newton's constant for the theory (|2.ip . 

Since Newton's constant G is dimensionful, we need 
to define the system of units used in measurements of 
G, as our final result will depend on the system of units 
chosen. Here we are not concerned with changes of units 
in the usual sense where the ratios between the old and 
new standards of mass, length and time are constants, 
independent of space and time. Rather, we are concerned 
with changes in the operational procedures for how units 
are defined, for which the ratios between the old and new 
standards can vary with space and time, as discussed 
by Dicke [i^. Changes of units of this type encompass 



Note that the critical value of velocity dispersion cr = /3 is a 
factor of \/3 larger than indicated by the approximate formula 
jBgl which was used in Ref. llSll. 
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conformal transformations of the metric, but can be more 
general. The necessity of specifying a system of units 
before discussing the space-time variation of dimensionful 
constants of nature is discussed in detail in Duff [43 |. 

We choose to use systems of units which are deter- 
mined completely by physics in the visible sector (which 
we label by j = 6 for baryonic), and are not determined 
by gravitational physics or by physics of the dark matter 
sector. For example, this is true of SI units. Alterna- 
tively one could take one's standards of mass, length and 
time to be the mass, size, and light travel time across a 
Hydrogen atom. Any system of units of this type will 
give rise to the formula (jCip below for G. 

We note that alternative choices arc in principle pos- 
sible. For example, one could have a small fiducial black 
hole, and one could transport it adiabatically around the 
Universe to act clS cL standard. Then one could define 
the standard of mass to be the mass of the black hole, 
while defining the standards of length and time in terms 
of Hydrogen atoms as above. In this system of units G 
would vary in space and time, but in a manner different to 
that described by the formula (jCip below, and Planck's 
constant h would also vary in space and time. Another 
possibility would be to define the units of mass in terms 
of the mass of the dark matter particle; this would also 
yield a different formula for G whenever the dark matter 
coupling function q:c{4>) differs from that of the baryons 
abi4>)- 

We next discuss the quantities on which the formula for 
G depends. First, for the theory (|2.ip with several differ- 
ent matter sectors, Newton's constant becomes a matrix 
whose i,j element governs the strength of the gravita- 
tional interaction between sector i and sector j. Second, 
G will depend on the scalar field (j), and through this de- 
pendence become a function of space and time. Here 4> 
should be thought of as a background value of the scalar 
field; gravitational interactions that act as perturbations 
to this background, for which the perturbation to the 
scalar field can be treated linearly, have a strength de- 
scribed by the formula (|C1[) . Finally, since the scalar con- 
tribution to the gravitational force will in general have 
a Yukawa profile, G will also depend on a spatial wave 
vector k, which is measured in the units discussed above. 

Our formula for Newton's constant is 



"'off (0.P) 



(CI) 



appeared before in the literature. The overall prefactor 
of e^"*" is well known, and the formula for the case of one 
matter sector has been previously derived in the con- 
text of cosmological perturbation theory [HI. The 1 in 
the square brackets describes the exchange of a graviton, 
and the second term in the square brackets describes the 
exchange of a scalar quantum, which couples to particles 
in the sector i with an amplitude proportional to a[{(f)). 
The scalar coupling term vanishes in the long wavelength 
limit fc ^ if the scalar field has a finite effective mass 

TOcff. 

We now turn to the derivation of the formula (jCip . For 
the derivation it will be convenient to use a different sys- 
tem of units, which are defined as follows. The standards 
of mass, length and time are defined to be the mass, size 
and light travel time across a small fiducial black hole, 
which is transported adiabatically around the Universe 
to act as a reference. Equivalent units^^ can be defined 
by demanding that the speed of light c and Planck's con- 
stant h be unity, and that lengths and times are measured 
using the Einstein frame metric gab- At the end of the 
derivation we will transform back to the visible-sector- 
based units. 

Starting with the action (|2.ip , we perform the following 
steps. We specialize the matter action in the jth sector 
to be that of a collection of point particles with masses 
nia] and worldlines a:;Qj'"(A): 



'gab 



(C3) 

We specialize to the Newtonian limit, so that the 
Einstein-frame metric gab is of the form — (1 + 2^)dt^ + 
(1 — 2^)dx^. We write the scalar field as 4> + 5(j), and 
expand to quadratic order in Scf) and <&. We assume that 
there is a background mass density p, so that the effec- 
tive mass of the scalar field is given by Eq. (jC2p . Finally, 
we compute from the action the energy E for a static 
configuration, up to an additive constant. This yields 



E = 



(C4) 



Here on the right hand side, G is a constant, related to 
by G = l/(87rmp), and the effective mass moff(0, p) 



is defined by 



''off 



(0,p) 



g"Kff(0,p) 



(C2) 



where the effective potential is given by Eq. (|2.1ip . In 
the adiabatic regime this effective mass reduces to the 
mass p.Sp . Several different aspects of this formula have 



where 



Pi 



(C5) 



^® By equivalent units we mean units that differ only by multipli- 
cation by constants. 
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Now integrating out $ and 5(f> yields 



E 



'j(0) 



4m2 



i.j Ql,a2,--- /3l,/32,■■■ 
^2„//^^„//^^„-mottI•^^ 



'a/3 



l + 2m2a;(0)aj^(0)e- 



(C6) 



where rJJ^ 



Ix«-x«| 



Now we transform back to the visible sector or Jordan 





frame units. The energy E, distances fj^^ and masses 
rha] in those units are given by 



E = e~"'''^'^^E, 



~1J 



(C7a) 
(C7b) 

(C7c) 



where the factor in square brackets on the right hand 
side of Eq. ((U7c)) is the mass as measured in the Einstein- 
frame units. Substituting these unit transformations into 
the energy expression (|C6[) and transforming from posi- 
tion space to momentum space yields an expression for 
the energy from which we can finally read off the effective 
Newton's constant (jCip . 
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